On the apparent lack of power in the CMB anisotropy at large angular scales 
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We study the apparent lack of power on large angular scales in the WMAP data. We confirm that 
although there is no apparent lack of power at large angular scales for the full-sky maps, the lowest 
multipoles of the WMAP data happen to have the magnitudes and orientations, with respect to the 
Galactic plane, that are needed to make the large scale power in cut-sky maps surprisingly small. 
Our analysis shows that most of the large scale power of the observed CMB anisotropy maps comes 
from two regions around the Galactic plane (~ 9% of the sky). One of them is a cold spot within 
~ 40° of the Galactic center and the other one is a hot spot in the vicinity of the Gum Nebula. If 
the current full-sky map is correct, there is no clear deficit of power at large angular scales and the 
alignment of the I = 2 and I = 3 multipoles remains the primary intriguing feature in the full-sky 
maps. If the full-sky map is incorrect and a cut is required, then the apparent lack of power remains 
mysterious. Future missions such as Planck, with a wider frequency range and greater sensitivity, 
will permit a better modeling of the Galaxy and will shed further light on this issue. 
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Over the last few years, detailed measurements of 
the anisotropies in the Cosmic Microwave Background 
radiation (CMB) have provided us a wealth of infor- 
mation. Recently the Wilkinson Microwave Anisotropy 
Probe (WMAP) [H, [2, |3, |i| mission released maps of 3 
years of data of the full sky in five different frequency 
bands. While the success of the standard inflationary cos- 
mology in explaining these observations combined with 
other astronomical data is remarkable, there are appar- 
ent discrepancies between prediction and observations on 
large angular scales. One of these discrepancies is the 
lack of correlated signal on large angular scales in cut-sky 
maps of WMAP data. It seems that masking the sky with 
a Galactic mask reduces the correlation on large scales. 
This feature was also clearly seen by COBE [5|]. Figure 
[T] shows the measured two-point correlation in full- and 
cut-sky maps together with the prediction of the best 
flt LCDM model. The flatness of the two-point corre- 
lation of the WMAP data on large angular scales was 
flrst quantified by the a posteriori statistic of for the 
first-year data and in more detail by for the 3-year 
data. It was reported that the large-angle correlations, 
for 9 > 60°, are unusually weak at the 99.85% confidence 
level and it is more siginificant in three year data than 
in the first-year data [7| . 

The two-point correlation function, C{6), of a CMB 
anisotropy map is the Legendre transform of its power 
spectrum C'l, 

C{0)^ —y{2l + l)CiPi{cos9). (1) 
1=1 

On large angular scales, C(Q) is sensitive to the low order 
multipoles and low-Z features in the angular power spec- 
trum are projected onto large-angle patterns in the two- 
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FIG. 1: The two-point correlation functions for WMAP data. 
Cut-sky maps have surprisingly low power on angular scales 
larger than 60°. This is consistently seen in the ILC map, 
TOH map and WMAP foreground reduced maps for Q, V and 
W frequency bands masked with kpO intensity mask. Two- 
point correlation function for full-sky ILC and TOH maps are 
also shown. Unlike masked maps, they are not flat on large 
angular scales. 



point correlation. But the disadvantage of using C(6) 
is that it is highly correlated at different angular scales. 
Therefore a covariance analysis is suggested for reliable 
statistical studies of the two-point correlation [13] ■ 

To quantify the lack of power on large scales, [6| sug- 
gested the sum of the square of the two-point correlation 

1 /2 

/_i C'^{9)dco%B. That is a simple and useful statistic 
but since C(&) is correlated on different angular scales, it 
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is better to use a covariance weighted sum of the square 
of the two-point correlations defined by the following four 
point statistic 



A{x) 



C{9)F^^{9, e')C{e')dcoa Od cos 6' , (2) 



in which F{9, 6') is the covariance between different bins 
in C{9) defined as 



F{e,e') = {C{e)c{e')). 



(3) 



where (• • •) represents the ensemble average and C{0) = 
C{e) - {C{e)). Note that A{l/2) in the limit of uncor- 
related C{9), will be the same as the S'i/2-statistic origi- 
nally proposed by Q . 

We use the following maps and masks for our analysis: 

• The full-sky WMAP Internal Linear Combination 
(ILC) map released by the WMAP team. 

• The three-year foreground-cleaned map of Tcgmark 
et al. 13 (TOH map) [25]. 

• WMAP foreground reduced maps for Q, V and W 
frequency bands, co-added from each differencing 
assembly using noise weighting (see [l|). 

• KpO intensity mask. This mask excludes the 
Galactic region (in which the K-band intensity is 
high) and also 0.6° around known point sources 
(^^23. 46% of the pixels from the maps) at r9 resolu- 
tion corresponding to HEALPix resolution Ngide = 
512. WMAP maps and the mask are available on 
LAMBDA [2^ as part of the three year data re- 
lease. 

All maps and masks are at r9 resolution from the three 
year data of WMAP. The two-point correlations are com- 
puted using eqnlT]and the covariance matrix F{9,9'), is 
computed from simulations of the best fit LCDM model 
for each mask. Computing A{x) for the above maps 
shows that A{x) for cut-sky maps is much smaller than 
for full-sky maps. This is expected because of the small 
values of C{9) on large angular scales in masked maps. 

To assign a statistical significance to the smallness 
of A{x) computed from the data, we perform compar- 
isons against 10,000 random simulations of statistically 
isotropic and Gaussian CMB anisotropy maps based on 
the best fit LCDM model to the WMAP data. Since 
A{x) is positive definite, the statistical significance of 
A{x) at each point is defined as the fraction of the 
simulations that have a smaller A{x) than the data, 
P{Asim{x) < Adata{x))- The results are shown in Ta- 
ble |l] in which we have made an a posteriori choice for 
a; to be X ~ 0.53 {9 — 58°). This is the point at which 
A{x) of the masked ILC map reaches its minimum. We 
will use A{l/2) as a short name for A{x = 0.53). 



60 80 100 120 140 160 180 




FIG. 2: Apparent lack of power on large angular scales in 
cut-sky ILC (top) may be more than just a low quadrupole. 
Removing the quadrupole does not change the shape of the 
two-point correlation much (middle). But removing Z = 2, 3 
together has a bigger effect (bottom). 



The full-sky maps are not anomalous; 8% of 10,000 
simulations have an A(x) smaller than that of the full- 
sky WMAP maps (see for a different statistic but 
similar conclusion). The odd shape of the full-sky cor- 
relation, and in particular the anticorrelation at 9 — n, 
is due to the fact that even multipoles are weaker than 
the odd multipoles at small I {I < 20). The suppression 
of even multipoles was studied by [lj| for the first-year 
maps of WMAP and it was shown that the sawtooth pat- 
tern of the angular power spectrum at small I was not 
statistically significant although visually striking. 

As opposed to the full-sky maps, masked maps have 
a small probability: only 0.73% of simulations had an 
A{x) smaller than that of ILC map masked with the kpO 
intensity mask. Q and W frequency band maps have the 
same probability and for the V band it is 0.69%. For 
the TOH map this number is 0.95%. This quantifies the 
smallness of the two-point correlation on large scales of 
cut-sky maps (see Table 11). These results qualitatively 
agree with the results of 7|. 

Is the low power in the cut-sky maps due to the 
low quadrupole? Removing the quadrupole changes the 
shape of the two-point correlation to some extent. But 
removing I = 2 and I = 3 multipoles together seems to 
have a bigger effect on the shape of the C(9) (Fig. [2]). 
Interestingly, the contribution from the 1 — 2,3 multi- 
poles to the two-point correlation on large scales is al- 
most equal to the contribution of the rest of multipoles 
but with a negative sign. That is 



{2l+l)CiPiicosi 

1=2,3 



\e>60° 



l>3 



(2l+l)CiPi(cost 



e>60°- 



(4) 
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TABLE I: The significance of the smallness of C{9) on large 
angular scales for different sets of maps and masks discussed 
in this paper. The third column is the fraction of 10,000 
simulations that has an A{x) smaller than the A{x) derived 
from the 3-year data of WMAP. This shows that a) cut-sky 
maps are discrepant on large scales while full-sky maps are 
not, b) masking the Galactic center and the region around 
the Gum Nebula in full sky CMB maps (8.9% of the sky), 
maskl, has almost the same effect as the kpO mask (23.5% of 
the sky), while applying the complimentary mask (the rest 
of the kpO mask with these two spots excluded, see Fig. |3]) 
does not suppress the large scale power much, and c) the 
quadrupole and the octupole happen to have the magnitude 
and orientation needed for making the two-point correlation 
in cut-sky maps very small on large angular scales. Excluding 
them from the maps or changing their orientation brings the 
large scale power back. 
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(kpO) 


6.50% 


ILC 




(maskl) 


0.78% 


ILC 




(kpO-maskl) 


12% 



"■A{x) is calculated at x 
cut-sky maps). 



0.53 (the minimum of the A{x) for 



To quantify this effect, we compute the A{x) for 
quadrupole and octupole removed maps. The simula- 
tions are also made from the best fit LCDM power spec- 
trum with C2 = C3 = 0, and 40% of simulations have 
an A{l/2) smaller than that of the cut-sky ILC map (see 
Table |T]) , which confirms the effect of the quadrupole and 
octupole on making the large scale power small, and indi- 
cates that the magnitudes of the low multipoles happen 
to be such that they cancel the contribution from the 
rest of the multipoles to the correlation function on large 
scales. 

The fact that smallness of two-point correlation on 
large scales happens in the cut-sky maps suggests that 
the orientation of the low multipoles relative to the 
Galactic plane may also play a role. To examine that, we 
keep the magnitude of the quadrupole fixed and change 
its orientation. This can be done by permuting the har- 
monic coefficients, a2m, of the full-sky ILC map in the 
following way 01 : 



~> ^20 7 0,21 



0-22 



0.21- 



A CMB anisotropy map constructed with this new 
quadrupole will have the same statistical properties as 
the ILC map (including the same power spectrum), but 



the orientation of the quadrupole in this map is different 
from that of the ILC map. We apply the kpO mask to this 
map and compute the two-point correlation, which ap- 
pears to be normal. The A(l/2) of this map is improved: 
the fraction of simulations that have smaller A{l/2) than 
this map increases from 0.73% to 6.5%. This proves that 
not only the magnitude of the low multipoles but also 
their orientation is responsible for making the large scale 
power of the cut-sky maps negligible. 

The fact that the large scale power of the WMAP data 
comes from the galactic plane is strange enough, but we 
can hunt it down further to see what parts of the Galac- 
tic plane contain most of the power. To do so we break 
the kpO mask into different sections to make new masks. 
A{\/2) is then computed for the ILC map masked with 
each of these new masks and compared to 10,000 sim- 
ulations of CMB anisotropy sky masked with the same 
masks. We see that most of the large scale power comes 
from two spots, a cold spot at —41° < 9 < 37° and a 
hot spot at 244° < 9 < 276°, both regions within the 
kpO intensity mask. We make a new mask to cover those 
two spots (maskl). This mask excludes 8.9% of the pix- 
els from the maps. C{9) computed from the ILC map 
masked with maskl is flat and close to zero on large an- 
gular scales and has a shape similar to the C{9) derived 
from the ILC map masked with kpO. Also A{l/2) is small 
for this map and P{Asim < ^masfei) = 0.78%. On the 
other hand the complementary mask to maskl, which is 
the rest of the kpO mask, has a small effect on the shape 
of C{9) on large angular scales and C{9) for the ILC map 
masked with this complementary mask is close to C{9) of 
the full-sky ILC map on large scales. Also A{l/2) is not 
small for the ILC map masked with this complementary 
mask and P{Asim < Akpo-maski) = 12% (see Table 
This mask has a small effect on the octupole but changes 
the orientation of the quadrupole. This is in agreement 
with the observation of [l2| that the quadrupole (unlike 
the octupole) has a non-negligible dependence on the ap- 
plied mask and foreground correction. 

Our analysis shows that most of the large scale power 
of CMB anisotropy maps comes from two regions around 
the Galactic plane. One of them is a cold spot within 
~ 40° of the Galactic center and the other one is a hot 
spot in the vicinity of the Gum Nebula. Although these 
are the potential suspects for foreground contamination 
Tsl 14 1, the discrepancy arises when one masks them out. 



It might be possible that the actual power of the CMB 
anisotropy is intrinsically low on large angular scales. 
There are a few theoretical models that predict suppres- 
sion of power on large scales. Spatially compact spaces 
are examples of these models. Although compact spaces 
with fine tuned matter (and total) densities can explain 
the suppressed large scale power of the WMAP data (see 
e.g. la, llTf). they fail to explain quadrupole-octupole 
alignment [16[, are not found by the circles- in-the-sky 
search |18| and are inconsistent with the bipolar power 
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FIG. 3: Masking these two spots (~ 9% of the sky) in the 
ILC map makes the two-point correlation of the resultant map 
smaU on angular scales larger than 60 degrees. 



spectrum analysis of the WMAP data [13, [3, [l^l ■ Other 
possibihties such as fine tuned Bianchi models are already 
ruled out because their model parameters are not con- 
sistent with the observed cosmological parameters and 
particularly because their total density is too low (see 
2l| for a complete review on these models). Other pos- 
sibilities such as homogeneous local dust-filled voids that 
can suppress the fluctuations due to the linear ISW effect 
have recently been studied [2^. 

At present more data are needed to determine the na- 
ture of the regions that contain most of the large scale 
power. If the current foreground modeling is correct then 
the only remaining peculiarity is the I = 2 and / — 3 
alignment 2j|. If the current foreground modeling is 
incorrect and there is no CMB power on large angular 
scales, then we will have to find a better cosmological 
model or live with a ~ 1% (a posteriori) chance for our 
current preferred model. 
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map 



be 



found 



here: 



http: / /space. mit.edu/home/tegmark/wmap. html 
http://lambda.gsfc.nasa.gov/ 

This is not a general way of rotating the aim coefficients, 
but is an example of a possible way of changing the orien- 
tation of a multipole by preserving the Ci. A general way 
of doing so is by the means of Wigner rotation matrices 
(see e.g. JM])- 



